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Fast simulation of the whole-sky CMB map in the presence 
of primordial magnetic field 
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ABSTRACT 

We present a novel method for generation of sets of the Cosmic Microwave Background 
(CMB) anisotropy maps, which reproduces the At ~ 2 correlations associated with 
Alfven turbulence. The method is based on the non-linear transformation of the CMB 
maps, which is obtained from the Monte Carlo simulation of the statistically isotropic 
Gaussian signal. Our method is computationally fast and efficient. 

We have applied two estimators (the cross-correlation estimator in multipole do- 
main for €-f 1, TO and l — l,m modes and circular phase moments) to test the statistical 
properties of derived maps. Both of these statistics confirm the effectiveness of our gen- 
eration method. We believe that our method can be useful for fast generation of the 
non-Gaussian maps in the presence of the primordial magnetic field, and be a valuable 
tool for the non-Gaussianity investigation of the CMB in the framework of the future 
PLANCK data analysis. 
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1 INTRODUCTION 

The primordial magnetic field (PMF) could be one of the 
most impressive non-Gaussian relic of the epoch of infla- 
tion (see Branderburg, Enqviest and Olesen (1996), Dur- 
rer, Kahniashvili and Yates (1998), Mack, Kahniashvili 
and Kosowsky (2002), Subramanian, Seshadri and Barrow 
(2003), Chen et al. (2004), Naselsky et al. (2004), Giovan- 
nini (2007), and etc.) which can be tested by observation of 
CMB anisotropy and polarization. 

The presence of the PMF with present strength ~ 1 — 10 
nC and coherent scale comparable or above the present day 
horizon can induce and support vorticity or Alfven turbu- 
lence before and during the epoch of hydrogen recombina- 
tion and via interaction with the CMB photons can produce 
non-adiabatic tail of the CMB anisotropy and polarization. 
The existence of a large-scale coherence of the PMF is a 
potential source of the preferred angular directions in the 
CMB map, which manifest themself as a sources of statistical 
anisotropy and non-Gaussianity of the CMB. Complimen- 
tary to the models of non-trivial topology such as Bianchi 
Vllh, cosmological strings and etc., the PMF can signif- 
icantly contribute to the low multipole part of the CMB 
power spectrum, causing some morphological peculiarities 
of the CMB sky. 

As shown in Chen et al, 2004, the Alfven turbulence 
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arising from primordial magnetic field (statistically homo- 
geneous and isotropic at very large scales) , should have non- 
Gaussian properties, because of quadratic dependence of the 
vorticity amplitude on the magnetic strength B. In the refer- 
ence coordinate system of z||B, magnetic field induces corre- 
lation between the a;_i,m and a;+i,m harmonic coefficients 
of the CMB temperature anisotropy ^ . Hereafter, our choice 
of the reference coordinate system (z||B) is implicit through 
this paper. 

Since the PMF has drawn very serious attention after 
the WMAP experiment, the investigation of possible obser- 
vational traces from the PMF has been included in the data 
analysis agenda for the PLANCK surveyor mission. In the 
framework of the non-Gaussianity investigation of the CMB 
observation by the PLANCK we need to generate about 
10'^ — 10^ realizations of the CMB anisotropy maps in the 
presence of the PMF in order to test their statistical prop- 
erties by various methods. Keeping "the open window" for 
other sophisticated methods, in this paper we propose a fast 
generation method of the CMB PMF map, based on non- 
linear transformation of the statistically homogeneous and 
isotropic Gaussian Random Field into non- Gaussian and 
anisotropic one. 

Let us outline the basic requirements for any methods 
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of modeling the CMB PMF. Firstly, the method should cor- 
rectly reproduce the statistical properties of the CMB PMF 
signal, namely, the power spectrum and the correlations be- 
tween a;_i,m, and a;+i_m, harmonic coefficients. Secondly, the 
method should be computationally simple and fast prefably 
at the level of Random Gaussian Field (RGF) simulation. 
Thirdly, it should possess good scalability with increasing 
angular resolution of the maps. For instance, the method 
should cover the multipole range up to Z ~ 4 — 6 x 10^ for the 
PLANCK mission. The outline of this paper is as follows. 
In Section 2 we briefly discuss the basic properties of the 
CMB PMF. Section 3 is devoted to description of the non- 
linear generator (hereafter PMF generator) and discusses 
the properties of derived signal. In Section 4 we present the 
analysis of correlation between different multipole moments 
and show that our PMF generator reproduce correct theo- 
retical properties of the CMB PMF signal. In Section 5 we 
make the phase analysis on the generated CMB PMF maps 
in order to show that morphology of generated maps is cor- 
rectly reproduced by our PMF generator. In conclusion we 
summarize the results of investigations, and mark problems 
relavant to the future PLANCK data analysis. 



Moreover, since the statistical isotropy of the CMB has been 
broken by the presence of the PMF, the power spectrum 
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now reflect directly this statistical anisotropy through the 
dependence of Ce on m. 

In the paper Durrer, Kahniashvili and Yates (1998) 
(hereafter DKY) it was shown that the power spectrum 
Ciijn) and auto-correlator Di(m) are given by 
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2 STATISTICS OF THE CMB ANISOTROPY 
GENERATED BY THE ALFVEN 
TURBULENCE 

One of the method to characterize the temperature 
anisotropy in the direction n = (6', ip) is decomposition of 
the signal AT(n) in spherical spherical harmonics: 
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where \airn\ and (j)im arc the modulus and phase of the ex- 
pansion coefficients respectively. 

Homogeneous and isotropic CMB Gaussian random 
fields (GRFs) possess I, m modes whose real and imaginary 
parts are Gaussian and mutually independent. The statisti- 
cal properties of the RGF are then completely specified by 
its angular power spectrum C1"^^ , 
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In other words, from the Central Limit Theorem their phases 
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are randomly and uniformly distributed in the range [0, 2tt], 
and the amplitude \aim\ follows Ryley distribution. We here- 
after denote the pure Gaussian CMB signal by the super- 
script " G" . The bracket in Eq. 2 denotes average over an en- 
semble of realizations. The expectational value of the power 
spectrum over an ensemble corresponds to theoretical pre- 
dictions. For the magnetized Universe, the helical part of 
vorticity does not contribute to the power spectrum, but in- 
duces off-diagonal correlations between £—l,m and £ + l,m 
harmonic coefficients: 
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where n is the power index of the Alfven turbulence power 
spectrum and Eq. 6 is valid in the range —7 < n < —1. 
The simplest way to detect the presence of a homogeneous 
magnetic field in the Universe is based on the arithmetic 
means over m of the two spectra, proposed in DKY 
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DKY find that 
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where Aq is a dimensionless normalization constant, tdec/to 
is the decoupling time to the present time ratio, and fco is the 
damping wave number. Though it is not shown explicitly in 
Eq. 8, the values of Di{m) may also be negative. As one can 
see from Eq.(8), the statistical isotropy of the CMB signal 
from the PMF is now restored, since the power spectrum Ct 
does not depend on the m. Thus, if we model the statistical 
properties of the CMB PMF by using Ce and De powers as 
(Chen et al, 2004), we are actually dealing with statistically 
isotropic, but non-Gaussian random field. This model is very 
useful for the situation, where the coherence scale of the 
PMF is much smaller than the present-day horizon, and the 
patterns of magnetic field are randomly oriented in space 
just as magnetic domains in ferromagnetic materials. 



Using the limiting values for Ao and fco in Eq. 8, the 
foUowings Ci are obtained (Chen et al. 2004): 
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where B is the strength of magnetic field in the unit of 
Gauss. These two cases correspond to a Harrison-Peebles- 
Yu-Zel'dovich scale-invariant spectrum and a possible in- 
flation model of primordial vorticity field respectively. The 
WMAP data impose the foUwing upper limit on the strength 
of the primordial magnetic field at 3 ct confidence level (Chen 
et al. 2004): 



B < 15 X lO^^G (n = -5), 
B<1.7xlO"^G (n = -7). 
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3 NON-LINEAR GENERATOR OF THE CMB 
PMF. 

In this section we discuss the simplest way to model the 
properties of the CMB PMF, which uses RGF CMB signals 
obtained from the Monte Carlo method. Let's start from the 
model, which can reproduce the average power spectra Ce 
of the CMB PMF signal. For this, we generalize the method, 
proposed in Naselsky et al. (2006) and define a CMB map 
with the following harmonic coefficients:. 
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Note that in Eq. (11) the a-parameter is a function of the 
power index n only, which can be easily shown by Eq. 8 
and 18. We draw ae^rn in Eq. (11) randonly from Gaussian 
distribution of the variance Coif)- From Eq. (11) one can 
get 
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Since ae,m. corresponds to RGF, average over realizations 
and m gives us 
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Thus, the power of the RGF Co {£) is equivalent to the power 
of the CMB PMF, up to the factor l + a'^. Let's draw our at- 
tention to the following combination of the ce,m coefficients: 

Q-l,mQ+l,m + Q-l,mQ+l,m = 
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Since the average over realizations gives us 
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Table 1. the value of a for various n 
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-4 ±2.36 ±(1.18 ±0.062) 

-5 ±2.26 ±(1.13 ±0.53) 

-6 ±2.21 ±(1.1 ±0.47) 

-7 ±2.17 ±(1.09 ±0.42) 



by averaging over m we finally obtain: 
De ~ aCoie.). (16) 

We introduce 7 to denote the Ci to D; ratio: 
Ci 
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Using Eq. 8, 13 and 16, we obtain the following equation for 
the Q-parameter: 
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and corresponding solutions: 
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In Table 1, we show the values of 7 and a for various n. 
The generation procedure is quite straightforward. We draw 
the set of ae^m from Gaussian distribution of the variance 
Ce/{i + a^)- Then, we obtain ce,m by applying the non- 
linear generator in Eq. 11 to the set of the RGF ae,m- We 
present illustrative cases in Fig. 1 and 2, which are ob- 
tained up to multipoles I ^ 500, using the upper limit val- 
ues of B in Eq. 10. The top panel in Fig. 1 shows the 
RGF CMB map drawn from Gaussian distribution of the 
power spectrum C{£) oc £~ . The middle panel shows the 
map we have obtained, using Eq. 11 with n — —5 and the 
maximal positive root a — 1.66 from Eq. (18). The bot- 
tom plot is the same as the middle one except for the sign 
of a = —1.66. These figures clearly show that our PMF- 
genorator significantly changes the morphology of the input 
RGF map, inducing significant Ai = 2 correlations. Fig. 2 
illustrates the properties of PMF-generator for the model 
with C{£) oc £~*, and n = —7. It should be noted that for 
this figure we used the same RGF CMB with that of Fig.l 
and simply rescaled the power by implementation of the fil- 
ter P{£,n = -7,n = -5) = C{£,n = -7)/C{£,n = -5) 
(see Novikov et al, 2001). It is not surprising that for the 
n = ~7 model the low muhipole tail of the CMB PMF im- 
ago is now dominant over high multipoles and the image has 
significant large angular scale modulation. From Fig. 1-2 one 
can see common feature of our PMF generation method. For 
a positive a the zone 9 = 7r/2 contains some stripes, while 
zones around polar cups retain much of morphology similar 
to the input map. For a negative a we have opposite ten- 
dency - the zones around the polar cups are modified most. 
However, to characterize these distortion quantitatively we 
need to use appropriate estimators of non-Gaussianity and 
statistical anisotropy. Below we use the simplest two estima- 
tors, which will directly reflect the coupling between A£ = 2 
multipoles. 
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Figure 1. The input map for the RGF CMB with n = —5 power 
spectrum (top). Middle - the CMB PMF for a = 1.66. Bottom- 
the CMB PMF for a = -1.66. 



Figure 2. The input map for the RGF CMB with n = —7 power 
spectrum (top). Middle - the CMB PMF for a = 1.51. Bottom- 
the CMB PMF for a = -1.51. 



4 STATISTICAL PROPERTIES OF THE CMB 
PMF 

4.1 Cross-correlation coefficient between 
multipole momentum. 

To characterize the coupling between I, m and £+2, m modes 
of the CMB PMF we propose to use the coefficient of cross- 
correlation, defined in multipole domain as 
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This coefficient K{l,A£) is directly related to the cross- 
correlation between any multipole moment i and £ + A£, 
and can be used for the analysis of A£ = 2 correlation of 
the CMB PMF. On the other hand, it should be noted that 
the K{1, Af)-estimator for uncorrelated signal may possess 
small but non-zero values, since the /("(Z, A^)-estimator is 
appUed to a single realization of the CMB PMF. However, 
we expect to find the cross-correlation of A£ = 2 unusu- 
ally high in comparison with correlations of A£ = 1,3,5.... 
Since the nominator and the denominator in Eq.(19) are 



similar to De and Ce when averaged over realizations, the 
values of K{1, At = 2)-estimator is expected to be ~ I/7. 
However, for each single realization we should have statis- 
tical deviation of K{l,A£ = 2) from the expectation value. 
Nevertheless, it is important to note that the K{1, A£ = 2)- 
coefficient is a good measure of A£ — 2 cross-correlations. 
In Fig. 3, significant A£ = 2 correlations of the simulated 
CMB PMF are clearly shown in contrast to relatively negli- 
gible A£ = 1 correlations. As shown in Fig. 3, the coefficient 
K(l, A£ = 2) of increasing £ is getting closer to Dt/Ct, while 
K(l,A£ = 1) — > 0. It might be also noticed that the over- 
all magnitude of K{1, A£ = 2) is slightly lower than Dt/Ci. 
It is attributed to approximations we made in derivation of 
Eq. 16. Thus, in terms of the cross-correlation coefficient our 
PMF-generator reproduces the statistical properties of the 
CMB PMF signal pretty well except for the low multipoles, 
where the sampling variance is large. 



4.2 The correlation of phases of CMB PMF 

The second statistic we would like to propose is based on 
measuring the Al = 2 phase correlation. The basic idea 
is to introduce some special functions of phases, which get 
the minimum contribution from uncorrelated Gaussian tail 
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Figure 3. The coefficient of cross-correlation K{1, A£). The thick 
solid lines correspond to the K(l,A£ = 2) with n = —5, a = 
1.66 (the upper curve) and a = —1.66 (the lower curve). The 
cross and plus signs correspond to K{l,A£ =1) with a = ±1.66 
respectively. The dashed lines correspond to ±Di/C'i. 
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Figure 4. C{e, At) for CMB PMF phases. The solid Unes corre- 
spond to C(£, A^ = 2) with n = —5, a = 1.66 (the upper curve) 
and a = —1.66 (the lower curve). The cross and plus signs corre- 
spond to C(l-, Al = 1) with a = 1.66 and a = —1.66. Two dashed 
lines represent the limit when (ji^+A^.m = 4't m (the top line), and 
'l>l+Ai,m = <i>e,m + 7f(the bottom line). 



and the maximum contribution from non-Gaussian tail of 
the phases. For this purpose the simplest trigonometric mo- 
ments statistics (Fisher, 1993) seem to be very useful (see 
Naselsky et al, 2005). Let's define the following trigonomet- 
ric moments: 



C{£, Ae) = — 2_^ COs{(j>e,,n - (f)l+Ae,m), 
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S{£, Ae) = — 2J sin{(j>e^,n - (pi+Ae.m), 



(20) 



where (j)t,m denotes the phase of a^,m- However, because of 
finite number of m modes, A£ — 2 cross-correlations may 
exist spontaneously especially in the low multipole range. 
Order-of-magnitude estimate on this effect is as follows. For 
pure Gaussian signal of uncorrelated phases (pi^m one can 
obtain 
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On the other hand, C(£,AV) estimator has an asymptotic 
form ~ £2 in the limit of complete phase correlation, while 
S(£) -^ 0. We have applied C(£,Al) estimator to the CMB 
PMF signal we had analyzed previously with a K(£, A£) 
statistic. The values of Q(£,AV) are shown in Fig. 4, which 
clearly shows that A£ = 2 correlation is sigrnficant in com- 
parison with A£ = 1 correlation just as the K(£^ A£) analy- 
sis. As it is seen from non-zero values of A£ — 1 correlation, 
there exists spontaneous correlation due to sample variance, 
which is biggest at the lowest multipole. 



5 MODIFICATION OF THE PMF ESTIMATOR 

5.1 flipping signs of correlations 

Our CMB PMF generator by Eq.(ll) is designed in such 
a way that the parameter a is a constant of a fixed sign. 
However, a parameters may not be a constant of a fixed sign. 
For instance, the sign of a parameter in Eq.(ll) may be a 
function of £, depending on the PMF model. The simplest 
way to extend Eq.(ll) is to assume that a{£) — ao(— 1) • In 
these models, the sign of Al — 2 correlation is alternating 
through multipoles. One can easily see that the signs do not 
need to follow the regular (deterministic) rules such as the 
one above, but rather stochastic ones. We show the CMB 
PMF generated under deterministic rules in Fig. 5, where 
we may see how the morphology of the maps can change. 
For even multipole £, K(£, A£ = 2) of these CMB PMF 
signals correspond to the upper solid line in Fig. 3, while 
for odd multipole £, it corresponds to the lower solid in Fig. 
3. One can see from Fig. 5 that a{£) of alternating sign 
leads to more homogeneous morphology of the maps, while 
preserving the Al — 2 correlations. 



5.2 "Brute force" magnetism 

As shown in Eq. 8, De is comparable to Ce when n > —1. 
Let's consider the generation of the CMB PMF for the model 
where the relative significance of A£ = 2 is highest (i.e. 
7 ~ 1). On the other hand, for 7 < 2, Eq. 18 fails to yield a 
real- valued a. Extending a to be a complex number requires 
Eq. 17 to be rewritten as follows: 



l-7Re[Q]4-|Q|^ =0. 



(22) 



However, it is not difficult to show that Eq. 22 is not satisfied 
by Q of any complex value either, if 7 < 2. Hence, the CMB 
PMF generator shown in Eq. 11 is not adequate for n > — 1, 
and we propose the following modified generator for n > — 1: 
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correlations resembles residual foregrounds, localized at the 
Galactic plane. Thus, our "brute force magnetism" method 
may be also used for modelling 1// noise or Galactic fore- 
grounds. 

Since this "brute force" technique impose full phase co- 
herence over AZ = 2 multipoles, the trigonometric moments 
of this "brute force" map correspond to the dash lines in 
Fig. 4. From Eq. 23, we may see that phases of CMB PMF 
satisfy the foUowings: 



^ — 2n,m — 



-I- TITV, lo = I ~ 21 



(24) 



where n is an arbitrary integer and [ J denotes the smallest 
integer larger than the argument. 

It should be noted that in the "brute force magnetism" 
the phases of the quadrupole component determine all the 
phases of even multipoles of |?n| ^ 2, while the phases of 
the octupole determine the phases of all odd multipoles of 
\m\ ^ 3. 



Figure 5. The PMF maps obtained with a{£) = ao(—l)^, ag 
1.66 (top) and qq = —1.66 (bottom). 




6 CONCLUSION 

We have presented the fast generation method of the CMB 
PMF, which is based on the non-linear transformation of the 
Monte Carlo simulated Gaussian map. Our method is com- 
putationally fast and efficient, and possesses good scalability 
with increase in angular resolution of simulated maps. The 
generated maps are non- Gaussian and satisfy the A£ = 2 
correlations of the CMB anisotropy generated by Alfven tur- 
bulence. 

We have tested the statistical properties of our CMB 
PMF maps by estimating the cross-correlation and circu- 
lar phase moments. As shown in this paper, both of these 
statistics have proved their effectiveness of our generation 
method. It turns out that our original generation method 
fails for Alfven turbulence of power index n > — 1, where 
De ~ Ci. For such Alfven turbulence, we have developed a 
so-called "brute force magnetism" method, whose generated 
maps possess full phase coherence over A£ = 2 multipoles. 
Our method can be easily extended to polarized signal and 
incomplete sky coverage (see Chiang and Naselsky, 2007), 
whose discussion will follow in separate publications. We 
believe our method is quite useful for generation of non- 
Gaussian maps associated with the primordial magnetic field 
and it will be a valuable tool for the non-Gaussianity study 
of the CMB in the framework of the future PLANCK data 
analysis. 



Figure 6. The map for positive Ai = 2 correlations (top) and 
negative correlation (bottom), obtained with Eq. 23. 



where ai^^n is drawn from Gaussian distribution of variance 
Ci and plus (minus) sign in front of the right hand side cor- 
responds to positive (negative) A£ = 2 correlations respec- 
tively. In Fig. 6 we show the CMB PMF maps generated 
by Eq. 23. The top map corresponds to positive Al — 2 
correlations and the bottom map corresponds to negative 
ones. As one can see from Fig. 6, the top map of positive 
AZ = 2 correlations resembles a set of longitudinal stripes 
and have resemblance to the properties of l//-noise widely 
discussed in the literature. The map of negative Al = 2 
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